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Q Abstract 

We consider multicommodity flow and cut problems in polymatroidal networks where there are sub- 
modular capacity constraints on the edges incident to a node. Polymatroidal networks were introduced by 
Lawler and Martel [24] and Hassin [ 18 1 in the single-commodity setting and are closely related to the sub- 
^ modular flow model of Edmonds and Giles 1 10|; the well-known maxflow-mincut theorem holds in this 

C/^ more general setting. Polymatroidal networks for the multicommodity case have not, as far as the authors 

are aware, been previously explored. Our work is primarily motivated by applications to information flow 
^ in wireless networks. We also consider the notion of undirected polymatroidal networks and observe that 

Q they provide a natural way to generalize flows and cuts in edge and node capacitated undkected networks. 

We establish poly-logarithmic flow-cut gap results in several scenarios that have been previously con- 
I sidered in the standard network flow models where capacities are on the edges or nodes ||25ll26l[T4l|23l[T3]| . 

^ Our results from a preliminary version have already found applications in wireless network information 

CNj flow ||20] l2n and we anticipate more in the future. On the technical side our key tools are the formula- 

tion and analysis of the dual of the flow relaxations via continuous extensions of submodular functions, 
in particular the Lovasz extension. For directed graphs we rely on a simple yet useful reduction from 
polymatroidal networks to standard networks. For undirected graphs we rely on the interplay between the 
Lovasz extension of a submodular function and line embeddings with low average distortion introduced 
by Matousek and Rabinovich [29 J ; this connection is inspked by, and generalizes, the work of Feige, Haji- 
aghayi and Lee ifTSll on node-capacitated multicommodity flows and cuts. The applicabiUty of embeddings 
to flow-cut gaps in polymatroidal networks is of independent mathematical interest. 
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1 Introduction 



Consider a communication network represented by a directed graph G = {V,E). In the so-called edge- 
capacitated scenario, each edge e has an associated capacity c(e) that limits the information flowing on it. We 
consider a more general network model called the polymatroidal network introduced by Lawler and Martel 
1241 and independently by Hassin [ 18 1. This model is closely related to the submodular flow model introduced 
by Edmonds and Giles lilOJ . Both models capture as special cases, single-commodity s-t flows in edge- 
capacitated directed networks, and polymatroid intersection, hence their importance. Moreover the models 
are known to be equivalent (see Chapter 60 in ||35]| . in particular Section 60.3b). The polymatroidal network 
flow model is more directly and intuitively related to standard network flows and one can easily generalize it 
to the multicommodity setting which is the focus in this paper. 

The polymatroidal network flow model differs from the standard network flow model in the following 
way. Consider a node u in a directed graph G and let 5q{v) be the set of edges in to v and 6q{v) be the set of 
edges out of v. In the standard model each edge (u, v) has a non-negative capacity c(n, v) that is independent 
of other edges. In the polymatroidal network for each node v there are two associated submodular functions 
(in fact polymatroid^ and p+ which impose joint capacity constraints on the edges in 6q{v) and 5q{v) 
respectively. That is, for any set of edges S C 5q{v), the total capacity available on the edges in S is 
constrained to be at most p^{S), similarly for 5q{v). Note that an edge [u, v) is influenced by /?+ and p~ . 
Lawler and Martel considered the problem of finding a maximum s-t flow in this model. The results in | l24l[T8l 
show that various important properties that hold for s-t flows in standard networks generalize to polymatroid 
networks; these include the classical maxflow-mincut theorem of Ford and Fulkerson (and Menger) and the 
existence of an integer valued maximum flow when capacities are integral. 

The original motivation for the Lawler-Martel model came from an application to a scheduling problem 
ll28l . More recently, there have been several applications of polymatroid network flows, (and submodular 
flows) and their generalizations such as linking systems Il36l . to information flow in wireless networks HI 12 
[38l[I21[33l|20l. A node in a wireless network communicates with several nodes over a broadcast medium and 
hence the channels interfere with each other; this imposes joint capacity constraints on the channels. Several 
interference scenarios of interest can be modeled by submodular functions. Most of the work on this topic 
so far has focused on the case of a single source. In this paper we consider multicommodity flows and cuts 
in polymatroidal networks where several source-sink pairs (si, ti), (s2, t2), • • • , (sfc, t^) share the capacity of 
the network. In the communications literature this is referred to as the multiple unicast setting. Our primary 
motivation is applications to (wireless) network information flow; see companion papers [20, 21 ] that build on 
results of this paper. Another motivation is to understand the extent to which techniques and results that were 
developed for multicommodity flows and cuts in standard networks generalize to polymatroidal networks. We 
note that polymatroidal networks allow for a common treatment of edge and node capacities; an advantage 
is that one can define cuts with respect to edge removals while the cost is based on nodes. As far as we are 
aware, multicommodity flows and cuts in polymatroidal networks have not been studied previously. 

Flow-cut gaps in polymatroidal networks: The main focus of this paper is understanding multicommodity 
flow-cut gaps in polymatroidal networks. In communication networks cuts can be used to information theo- 
retically upper bound achievable rates while flows allow one to develop lower bounds on achievable rates by 
combining a variety of routing and coding schemes. Flow-cut gaps are of therefore of much interest. Unlike 
the case of single-commodity flows where maximum flow is equal to minimum cut, it is well-known that even 
in standard edge-capacitated networks no tight min-max result holds when the number of source-sink pairs 
is three or more (two or more in case of directed graphs). See |[35l for some special cases where min-max 

'a set function / : 2^ ^ M over a finite ground set is submodular iff f{A) + f{B) > f{A n B) + f{A U B) for all 
A,B (ZN; equivalently f{A U {i}) - f{A) > f{B U {i}) - f{B) for all A C B and i ^ A. It is monotone if f{A) < f{B) for 
all A C B.ln this paper a polymatroid refers to a non-negative monotone submodular function with /(0) — 0. 
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results do hold. Flow-cut gap results have been extensively studied in theoretical computer science starting 
with the seminal work of Leighton and Rao ||25]| . The initial motivation was approximation algorithms for 
cut and separator problems that are NP-Hard. There has been much subsequent work with a tight bound of 
0(logA;) established for flow-cut gaps in undirected graphs in a variety of settings |[T4l l26l l6l [T3l . It has 
also been shown that strong lower bounds exist for flow-cut gaps in directed graphs; for instance the gap is 
0(min{A;, n^}) between the maximum concurrent flow and the sparsest cut ll34l l9l where 5 is a fixed con- 
stant. However poly-logarithmic upper bounds on the gaps are known for the case of symmetric demands in 
directed graphs ll23l[TT.l . Motivated by the above positive and negative results we focus on those cases where 
poly-logarithmic flow-cut gaps have been established. We show that several of these gap results extend to 
polymatroid networks. Our results and techniques lead to new approximation algorithms for cut problems in 
polymatroidal networks which could have future applications. However, in this paper we restrict our attention 
to quantifying flow-cut gaps. 

Bidirected and undirected polymatroidal networks: As we mentioned already, strong lower bounds exist 
on flow-cut gaps for directed networks. Positive results in the form of poly-logarithmic upper bounds on flow- 
cut gaps for standard networks hold when the demands are symmetric or when the supply graph is undirected. 
A natural model for wireless networks is the bidirected polymatroidal network. For two nodes u and v it is 
a reasonable approximation to assume that the channel from u to w is similar to that from v to u; hence one 
can assume that the underlying graph G is bidirected in that if the edge (n, v) is present then so is {v,u). 
Moreover, we assume that for any node v and S C 6~{v), p^{S) = pt{S') where S' C 6~^{v) is the set of 
edges that correspond to the reverse of the edges in S. Within a factor of 2 bidirected polymatroidal networks 
can be approximated by undirected polymairoidai networks: we have an undirected graph G and for each node 
V a single polymatroid p„ that constrains the capacity of the edges Sciv), the set of edges incident to v. The 
main advantage of undirected polymatroid networks is that we can use existing tools and ideas from metric 
embeddings to understand flow-cut gap results. Undirected polymatroidal networks have not been considered 
previously. We observe that they allow a natural way to capture both edge and node-capacitated flows in 
undirected graphs. To capture node-capacitated flows we set Pv{S) = 2c{v) for all 7^ C S{v) where c{v) 
is the capacity of We mention an advantage of using polymatroidal networks even when considering the 
special case of node-capacitated flows and cuts: one can define cuts with respect to edges even though the 
cost is on the nodes. This is in fact quite natural and simplifies certain aspects of the algorithms in [|13,1 . 

1.1 Overview of results and technical ideas 

We do a systematic study of flow-cut gaps in multicommodity polymatroidal networks, both directed and undi- 
rected. Let G = {V,E) be a polymatroidal network on n nodes with A; source-sink pairs {si,ti), . . . , {sk,tk)- 
We consider two flow problems and their corresponding cut problems: (i) maximum throughput flow and 
multicut (ii) maximum concurrent flow and sparsest cut. Our high-level results are summarized below. 

• For directed networks we show a reduction based on the dual that establishes a correspondence between 
flow-cut gaps in polymatroidal networks and the standard edge-capacitated networks. This allows us 
to obtain poly-logarithmic upper bounds for flow-cut gaps in directed polymatroidal networks with 
symmetric demands via results in |[23l [TTI for both throughput flow and concurrent flow. In particular 
we obtain an 0(min{log'^ k, log^ n log log n}) gap between the maximum concurrent flow and sparsest 
cut. The reduction is applicable only to directed graphs. 

• We show that line embeddings with low average distortion ll29l [32l lead to upper bounds on flow-cut 
gaps in polymatroidal networks — this connection is inspired by the work in |[T3l for node-capacitated 

^The factor of 2 is needed since a flow patli p through an internal node v uses two edges. On the other hand it is not needed for 
the sources and sinks. This technical issue is a minor inconvenience with undirected polymatroidal networks; we note that this also 
arises in treating node-capacitated multicommodity flows L13J . 
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flows. For undirected polymatroidal networks this leads to an optimal 0(log k) gap between maximum 
concurrent flow and sparsest cut. We also obtain an optimal 0(logA;) gap between throughput flow 
and multicut. These imply corresponding results for bidirected networks, which have already found 
applications 1201 |2T1 . As in ||T3l the embedding connection can be exploited to obtain improved ap- 
proximation algorithms for certain separator problems by exploiting graph structure |[22l or by using 
stronger relaxations via semi-definite programming (and associated embedding theorems) H |2l ; we 
defer these improvements to a later version. 

Most of the literature on multicommodity flow-cut gaps is based on analyzing the dual of the linear 
program for the flow which can be viewed as a fractional relaxation for the corresponding cut problem. The 
gap is established by showing the existence of an integral cut within some factor of the relaxation. For 
standard edge and node-capacitated network flows the dual Unear program has length variables on the edges 
which induce distances on the nodes. The situation is more involved in polymatroidal networks, in particular. 



the definition of the cost of a cut is some what complex and is discussed in more detail in Section 2.2 Our 
starting point is the use of the Lovasz extension of a submodular function L27,| to cleanly rewrite the dual of 
the flow linear programs. This simplifies the constraint structure of the dual at the expense of making the 
objective a convex function. However, we are able to exploit properties of the Lovasz extension in several 
ways to obtain our results. Our techniques give two new dual-based proofs of the maxflow-mincut theorem 
for single commodity polymatroid networks that was first established by Lawler and Martel algorithmically 
Il24l via an augmenting path based approach. We believe that the applicability of embedding based methods 
for polymatroidal networks is of independent mathematical interest. 

For the most part we ignore algorithmic issues in this paper although all the flow-cut gap results lead to 
efficient algorithms for finding approximate cuts. 



2 Multicommodity Flows and Cuts in Polymatroidal Networks 

We let G = {V, E) represent a graph whether directed or undirected. We use (li, v) for an ordered pair of 
nodes and uv to denote an unordered pair. In a directed graph G, for a given node v, 6q (v) and 6q (v) denote 
the set of incoming and outgoing edges at v. In undirected graphs we use 5c{v) to denote the set of edges 
incident to v. We omit the subscript G if it is clear from the context. In addition to the graph the input consists 
of a set of k source-sink pairs (si,ti), . . . , {sk,tk) that wish to communicate independently and share the 
network capacity. 

In a directed polymatroidal network each node v ^ V has two associated polymatroids p~ and />+ with 
ground sets as 6^ (v) and 5~^{v) respectively. These functions constrain the joint capacity on the edges incident 
to V as follows. If 5 C 5~ (ti) then (5) upper bounds the total capacity of the edges in S and similarly 
if 5 C 5~^{v) then p^{S) upper bounds the total capacity of the edges in S. We assume that the functions 
Pv{')^Pv{') V ^ V w& provided via value oracles. In undirected polymatroidal graphs we have a single 
function p^{-) at a node v that constrains the capacity of the edges incident to v. Continuous extensions of 
submodular functions, namely the Lovasz extension fTT\ and the convex closure, are important technical tools 
in interpreting and analyzing the duals of the linear programs for multicommodity flow in the polymatroid 



setting. We discuss these in Section 2.2 We first discuss the two flow problems of interest, namely maximum 



throughput flow and the maximum concurrent flow. 



2.1 Flows 

A multicommodity flow for a given collection of k source-sink pairs (si, ti), . . . , (s^, tk) consists of k sepa- 
rate single-commodity flows, one for each pair (sj, tj). The flow for the i'th commodity can either be viewed 
as an edge -based flow fi : E ^ or as a path-based flow fi : Vi ^ M_|_ where Vi is the set of all simple 
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paths between Sj and t j in G. We prefer the path-based flow since it is more convenient for treating directed 
and undirected graphs in a unified fashion, and also to write the linear programs for flows and cuts in a more 
intuitive fashion. However, it is easier to argue polynomial-time solvability of the linear programs via edge- 
based flows. Given path-based flows fi,i = 1, . . . , for the k source-sink pairs, the total flow on an edge e 
is defined as /(e) = Yh^i Y^pav^-.p^e f^iP)- The total flow for commodity iis Ri = YpeV, fiiP) where Ri 
is interpreted as the rate of commodity flow i. In directed polymatroidal networks, the flow is constrained to 
satisfy the following capacity constraints. 

^f{e)<p;;iS) yv\fS^6-iv) and 5^/(e)</>+(5) V^;V5 C 

The constraints in M?i^//recfe^/ polymatroidal networks are: J2e€S f(^) — Pvi'^) Vf VS* C 6{v). 

A rate tuple {Ri, ...,Rk) is said to be achievable if commodities 1, . . . ,k can be sent at rates Ri, . . . , Rk 
simultaneously between the corresponding source-sink pairs. For a given polymatroidal network and source- 
sink pairs the set of achievable rate tuples is easily seen from the above constraints to be a polyhedral set. We 
let P{G, T) denote this rate region where G is the network and T is the set of given source-sink pairs. In 
the maximum throughput multicommodity flow problem the goal is to maximize Ya=i ^« ^^^^ P{G,T). In 
the maximum concurrent multicommodity flow problem each source-sink pair has an associated demand Di 
and the goal is to maximize A such that the rate tuple {XDi, XD^) is achievable, that is the tuple belongs 
to P{G, T). It is easy to see that both these problems can be cast as linear programming problems. The 
path-formulation results in an exponential (in n the number of nodes of G) number of variables and we also 
have an exponential number of constraints due to the polymatroid constraints at each node. However, one can 
use an edge-based formulation and solve the linear programs in polynomial time via the ellipsoid method and 
polynomial-time algorithms for submodular function minimization. 

Network with symmetric demands: In directed polymatroidal networks we are primarily interested in sym- 
metric demands: node Si intends to communicate with tj and node ti intends to communicate with s,; at the 
same rate. Conceptually one can reduce this to the general setting by having two commodities (sj, ti) and 
{ti, Si) for a pair Siti and adding a constraint that ensures their rates are equal. To be technically consis- 
tent with previous work we do the following. We will assume that we are given k unordered source-sink 
pairs siti, . . . , s^tk- Now consider the 2k ordered pairs (si, ti), . . . , (s^, t^), {ti, si), . . . , (t^, s^). We are 
interested in achievable rate tuples of the form {Ri, . . . , Rk, R'l, ■ ■ ■ , R'j^) where R[ = Ri. In the maximum 
throughput setting we maximize Y'i=i{^i + ^D- Note that even though the rates for {si,ti) and {ti,Si) 
are the same, the flow paths along which they route can be different. In the maximum concurrent flow set- 
ting both {si, ti) and (tj, Si) have a common demand Di and we find the maximum A such that rate tuple 
{XDi, ...,XDk, XDi, XDk) is achievable for the pairs (si,ti), . . . , (sfc,ife), {h,si), {tk,Sk). 

2.2 Cuts 

The multicommodity flow problems have natural dual cut problems associated with them. Given a graph G = 
{V, E) and a set of edges F <^ E we say that the ordered node pair {s, t) is separated by F if there is no path 
from s to t in the graph G[E \F]. In directed graphs F may separate (s, t) but not (t, s). In undirected graphs 
we say that F separates the unordered node pair st if s and t are in different connected components of G[E\F]. 
For certain problems, especially in the information theoretic setting, it is of interest to consider restricted cuts 
induced by vertex bi-partitions, that is cuts of the form F = Sq{S) {F = 6g{S) in the undirected setting) 
for some S* C y. In this paper we mainly consider edge-cuts; Section [5. 1 . 1 [ discusses how vertex bi-partition 
cuts can be obtained from edge-cuts for the sparsest cut problem in undirected polymatroidal networks. 

In the standard network model the cost of a cut defined by a set of edges F is simply J2eeF '^(^) where 
c(e) is the cost of e (capacity in the primal flow network) . In polymatroid networks the cost of F is defined 
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in a more involved fashion. Each edge {u, v) in F is assigned to either u or v; we say that an assignment of 
edges to nodes g : F ^ V is valid if it satisfies this restriction. A valid assignment partitions F into sets 
I V G V} where g^^{v) (the preimage of v) is the set of edges in F assigned to v by g. For a given 
valid assignment g of F the cost of the cut i^g{F) is defined as 

u,{F) :=Y,{p;;i5-{v)ng~\v)) + pt{S+iv)ng-\v))). 

V 

In undirected graphs the cost for a given assignment is J2v Pv{g^^{v))- 

Given a set of edges F we define its cost to be the minimum over all possible valid assignments of F to 
nodes, the expression for the cost as above. We give a formal definition below. 

Definition 1. Cost of edge cut: Given a directed polymatroid network G = {V, E) and a set of edges F (1 E, 
its cost denoted by y{F) is 

mill y^{pz{5-{v)r^g-\v)) + pU5^{v)r^g-\v))) . (1) 

g:F-^V, g valid ^ — ' 



V 



In an undirected polymatroid network ^{F) is 



mm 



s>y, g valid ^—^ 



Pv{g-\v)). (2) 



V 



Lemma 1. The cut cost function is sub-additive, that is, ^{F U F') < v{F) + u{F')for all F, F' C E. 

Although not obvious, v can be evaluated in polynomial time via an algorithm to compute an s-t maximum 
flow problem in a polymatroid network. We do not, however, rely on it in this paper. 
We now define the two cuts problems of interest. 

Definition 2. Given a collection of source-sink pairs . . . , {sk,tk) in G = (y,E) and associated 

demand values Di, . . . , D^, and a set of edges F <^ E the demand separated by F, denoted by D{F), is 
^v(si u) separatedhy F ^i- ^ ^ multicut if all the given source-sink pairs are separated by F. The sparsity of 
F is defined as j)^- 

The above definitions extend naturally to undirected graphs. Given the above definitions two natural 
optimization problems that arise are the following. The first is to find a multicut of minimum cost for a given 
collection of source-sink pairs. The second is to find a cut of minimum sparsity. These problems are NP-Hard 
even in edge-capacitated undirected graphs and have been extensively studied from an approximation point 
of view ||25l[Il|26llliail. 

Lemma 2. Given a multicommodity polymatroidal network instance, the value of the maximum throughput 
flow is at most the cost of a minimum multicut. The value of the maximum concurrent flow is at most the 
minimum sparsity. 

A key question of interest is to quantify the relative gap between the flow and cut values. These gaps 
are relatively well-understood in standard networks and the main aim of this paper is to obtain results for 
polymatroid networks. 

Network witli symmetric demands: For a directed network with symmetric demands the notion of a "cut" 
has to be defined appropriately. We say that a set of edges F separates a pair Sjtj if it separates {si,ti) or 
{ti, Si). With this notion of separation, the definitions of multicut and sparsest cut extend naturally. A multicut 
is a set of edges F whose removal separates all the given pairs. Similarly for a set of edges F its sparsity is 
defined to v{F) / D{F) where D{F) is the total demand of pairs separated; note that if both {si,ti) and {ti, Si) 
are separated by F we count Di twice in D{F). This is to be consistent with the definition of flows given 
earlier. Lemma [2] extends to the symmetric demand case with the definition of flows given for symmetric 
demands in the previous section. 



5 



3 Relaxations for Cuts 



Lemma [2] gives a way to lower bound the value of multicuts and sparsest cuts via corresponding flow prob- 
lems. The flow problems can be cast as linear programs. The duals of these linear programs can be directly 
interpreted as linear programming relaxations for integer programming formulations for the cut problems. 
Here we take the approach of writing the formulation with a convex objective function and linear constraints; 
this simplifies and clarifies the constraints and aids in the analysis. For one of the cases we show the equiva- 
lence of the formulation with the dual of the corresponding flow Unear program. We first discuss continuous 
extensions of submodular functions. 

3.1 Continuous extensions of submodular functions 

Given a submodular set function p : 2^ ^ M on a finite ground set N it is useful to extend it to a function 
p' : [0, 1]^ —7- M defined over the cube in |A^| dimensions. That is we wish to assign a value for each 
X e [0, 1]^ such that p'(ls) = p{S) for all S C where I5 is the characteristic vector of the set S. 
For minimizing submodular functions a natural goal is to find an extension that is convex. We describe two 
extensions below. 

Convex closure: For a set function p : 2^ — )• M (not necessarily submodular) its convex closure is a 
function p : [0, 1]^ — )• M with p(x) defined as the optimum value of the following linear program: 

p(x) = mill ^ asp{S) 

SCN 

S.t. 

Y^as = 1 

5 

as = Xi \/i £ N 

S:ies 

as > ys. 

The function p is convex for any p. Moreover, when p is submodular, for any given x, the linear program above 
can be solved in polynomial time via submodular function minimization and hence p(x) can be computed in 
polynomial time (assuming a value oracle for p). It is known and not difficult to show that if p is a polymatroid 
(monotone and /(0) = 0) the value of the linear program does not change if we drop the constraint that 
E5«S = 1- 

Lovasz extension: For a set function p : 2^ — )• R (not necessarily submodular) its Lovasz extension ||27]| 
denoted by p : [0, 1]^ — > M is defined as follows: 

p(x)= [' p{^')de 

Jo 

where = {i \ Xi > 9}. This is not the standard way the Lovasz extension is stated but is entirely equivalent 
to it. The standard definition is the following. Given x let zi, . . . , i„ be a permutation of {1,2,..., n} such 
that X j^ > Xi^ > . . . > Xi^ > 0. For ease of notation define xq = 1 and Xn+i = 0. For 1 < j < n let 
Sj = {ii,i2, • • • , ij}- Then 

n 

p(x) = (1 - Xijp(0) + ^ixi^ - Xi^^^)p{Sj). 
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It is typical to assume that p{^) = and omit the first term in the right hand side of the preceding equation. 
Note that it is easy to evaluate />(x) given a value oracle for p. 
We state some well-known facts. 

Lemma 3. For a submodular set function p, p(x) = p{'K) for any x S [0, 1]^. Therefore the convex closure 
coincides with the Lovdsz extension and p{-) is convex. 

Proposition 1. For a monotone submodular function pandit < x' (coordinate-wise), /5(x) < /5(x'). 

The equivalence of p and p also implies that an optimum solution to the linear program defining /)(x) is 
obtained by a solution a where the support of a is a chain on N (a laminar family whose tree representation 
is a path). In fact we have the following. Given x G [0, 1]^ consider the ordering of the coordinates and 
the associated sets as in the definition of the /5(x). One can verify that as^ = Xi- — Xi^_^ for 1 < j < n, 
00 = (1 — Xi^), and as = for all other sets S is an optimum solution to the linear program that defines 
/5(x). We will use this fact later. 



3.2 Multicut 

We now consider the multicut problem. Recall that we wish to find a subset F (1 E such that F separates 
all the given source-sink pairs so as to minimize the cost ^{F). The only difference between the polymatroid 
networks and standard networks is in the definition of the cost. We first focus on expressing the constraint that 
F is a feasible set for separating the pairs. For each edge e we have a variable l{e) G [0, 1] in the relaxation 
that represents whether e is cut or not. For feasibility of the cut we have the condition that for any path p from 
Si to ti (that is p G 'Pj) at least one edge in p is cut; in the relaxation this corresponds to the constraint that 
X^esp ^(^) — ^- Other words dist£(si, ti) > 1 where disti{u, v) is the distance between u and v with edge 
lengths given by i{e) values. 

We now consider the cost of the cut. Note that z^(-F) is defined by valid assignments of F to the nodes, 
and submodular costs on the nodes. In the relaxation we model this as follows. For an edge e = (n, v) 
we have variables £(e, u) and i{e, v) which decide whether e is assigned to u or v. We have a constraint 
£{e, u) + £(e, v) = £{e) to model the fact that if e is cut then it has to be assigned to either u or v. Now 
consider a node v and the edges in 6^{v). The variables £{e, v), e G 6'^{v) in the integer case give the set 
of edges S C 5^{v) that are assigned to v and in that case we can use the function Py{S) to model the cost. 
However, in the fractional setting the variables lie in the real interval [0, 1] and here we use the extension 
approach to obtain a convex programming relaxation; we can rewrite the convex program as an equivalent 
linear program via the definition of p. Let be the vector consisting of the variables £{e, v), e G S~{v) 
and similarly d+ denote the vector of variables £{e,v), e G 6^{v). The relaxation for the directed case is 
formally described in Fig[T]in the box on the left. For the symmetric demands case the relaxation is similar but 
since we need to separate either (sj, U) or (ti, Si) the constraint dist£(si, tj) > 1 is replaced by the constraint 
disti{si,ti) + disti{ti, Si) > 1. 



minJ^(p-(d-) + p+(d+)) (3) 

V 

£{e,u) +£{e,v) = £{e) e = {u,v)£E 

disti{si,ti) > 1 1 <i < k 

£{e),£{e, u),£{e, v) > e = {u,v) e E. 



miny^p^(d^) 

V 


(4) 


£{e,u) +£{e,v) = £{e) 


e = uv £ E 


disti(si,ti) > 1 


l<i<k 


£{e),£ie,u),£ie,v) > 


e = uv £ E. 



Figure 1: Lovasz-extension based relaxations for multicut in directed and undirected polymatroidal networks 
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For the undirected case we let d^, denote the vector of variables i{e, v), e S 6{v) and the resulting relax- 
ation is shown on the right in Fig[T] 

One can replace in the above convex programming relaxations by the convex closure; further, one 
can use the definition of p^ via a linear program to convert the convex program into an equivalent linear 
program. The resulting linear program can be shown to be equivalent to the dual of the maximum throughput 
flow problem. See Section|A]for a formal proof. 



3.3 Sparsest cut 

Now we consider the sparsest cut problem. In the sparsest cut problem we need to decide which pairs to 
disconnect and then ensure that we pick edges whose removal separates the chosen pairs. Moreover we are 
interested in the ratio of the cost of the cut to the demand separated. We follow the known formulation in 
the edge-capacitated case with the main difference, again, being in the cost of the cut. There is a variable yi 
which determines whether pair i is separated or not. We again have the edge variables i{e),i{e, u),£{e, v) to 
indicate whether e = {u, v) is cut and whether e's cost is assigned to u or v. If pair i is to be separated to the 
extent of yi we ensure that dist£(sj, tj) > t/j. To express sparsity, which is defined as a ratio, we normalize the 
demand separated to be Fig [2] has a formal description on the left for the directed case. For the symmetric 
demands case we have essentially the same relaxation; the constraint J2i ^jdist^(sj, ti) = 1 is replaced by 
the constraint J2i Di{diste{si,ti) + dist£(tj, Sj)) = 1. 



V 




+ /5+(d+)) 






= £(e) e = {u,v)eE 


k 

^ A •dist^( 

1=1 

i{e),e{e,u),l 




= 1 




>0 e = {u,v)eE. 



min^p^,(dj,) 

V 




i{e,u) + i{e,v) =e{e) 


e = uv G E 


k 




^ A • dist£(si,ti) = 1 




i=l 




i{e),i{e,u),i{e,v) > 


e = {u, v) G E. 



Figure 2: Relaxations for sparsest cut in directed and undirected polymatroidal networks 



The relaxation for the undirected case is shown on the right in Fig |2] where d„ is the vector of variables 

i{e,v),e e 6{v). 



4 Flow- Cut Gaps in Directed Polymatroidal Networks 

In this section we consider flow-cut gaps in directed polymatroidal networks. We show via a reduction that 
these gaps can be related to corresponding gaps in directed edge-capacitated networks that have been well- 
studied. We note that this reduction is specific to directed graphs and does not apply to undirected polyma- 
troidal networks. The embedding based approach for the undirected case that we discuss in Section [5] is also 
applicable to directed graphs. 

The reduction is similar at a high-level for both gap questions of interest and is based on the relaxations 
for the two cut problems that we described in Section [3] We take a feasible fractional solution for relaxation 
of the cut problem in question and produce an instance of a cut problem in an edge-capacitated network and 
a feasible fractional solution to the corresponding cut problem. We also provide a correspondence between 

'We need to argue that this leads to a vahd relaxation for sparsest cut; it follows from the fact that the polymatroid functions in 
the objective function are normalized (pi,(0) = 0) and in this case we have /3„(tx) = t/5„(x) for any scalar t £ [0, 1]. 
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feasible integer solutions to the edge-capacitated network instance and the original problem such that the cost 
of the solution is preserved. These correspondences allow us to translate known gap results for the edge- 
capacitated networks to polymatroidal networks. 



4.1 Details of the reduction 

Let G = (y, i?) be a directed graph and let £ : — )• M+ be a length function on the edges. We let dist£(n, v) 
be the shortest path distance from u to f in G with edge lengths £. Moreover, for each edge (n, v) let £(e, u) 
and £(e, v) be two non- negative numbers such that £(e) = £(e, n)+£(e, v). For a node v let d+ be the vector of 
£(e, v) values for all edges e € 5"^ {v) and similarly is the vector of £(e, v) values for edges in 5' {v). In the 
polymatroidal setting the cost induced by the edge length variables is given by X]Dev(P (^»7) + P^i'^t))- 
Note that for multicut we have that dist£(si, tj) > 1 for each demand pair (sj, U) while in sparsest cut we 
are interested in the ratio of the cost to Ylii ' dist^(si, ti). We now describe the construction of a graph 
H = {Vh, Eh) where Vh = V ^ V' (that is the nodes of G are also in H) and an edge length function 
f : Eh — )• M+ such that dist^(ii, v) = dist£'(u, v) for all u,v ^ V; that is the distances between nodes in V 
are the same in G and H. We also create an edge-cost (or capacity in the primal sense) function c : Eh K+- 
The construction will also establish the correspondence of cuts in G and H and their costs. 

The graph H = (V ^ V, Eh) is constructed as follows. To aid the reader we first describe the idea of 
the construction at a high-level. Consider a node v ^ V and the in-coming edges 6~ (v) and out-going edges 
6'^{v). In H we have nodes of V and build an in-tree T~ and an out-tree T+ that are rooted at v. The leaves 
of T~ are the edges in 6^{v) the leaves of are the edges in 6~^{v). Note that an edge {u,v) will thus 
participate in T+ and T~. Now for the formal details. The nodes of H, denoted by Vh, consist of the nodes 
y of G and additional nodes V'. V has two types of nodes. First, for each edge e S -E there is a node 7e. 
Second, for each node u € F we create two sets of nodes A^^ (v) and A^+(t;) where \N^{v)\ = n~ = \5q{v)\ 
and I {v)\ = = \6q{v)\', thus one node for each edge in 6~{v)yj 5~ {v) ; these will be the internal nodes 
of the trees T~ and r+ respectively. For notational convenience we refer to the j'th node in N^{v) as v~ 
and similarly for the j'th node in N^{v). 

Now we describe the edge set Eh of the graph H, the edge length function i' : Eh — and the 
cost function c : Eh — M+- The edge set is essentially prescribed by specifying the trees and T+ 
for each v £ V. Consider the vector d~{v) of values £{e,v) for e G ^g^'^)- Recall the definition of the 
Lovasz extension /5^(d^). We order the edges in 6~{v) as ei,e2, . . . , e^- where i{ej,v) > £{ej-^-l,v) for 
1 < j < and then p~{d~) = X^j(^(ej, — i{ej+i,v))p~ (Sj) where Sj = {ei, . . . , ej}. We associate 
the node v~ with the set Sj. The edge set of T,7 is defined as follows. For ease of notation we let 

J ^ n„ +1 

represent the node v. We create a directed path t;^ —)•... —)• v~_ — )• f ~_ = v with edge lengths 

Tin Til) -\-\ 

i'{v^,V2) = ^(ei, w) - i{e2,v),i'{v2 ,v^) = i{e2,v) - i{e3,v), . . . ,i'{v~_,v) = £{e^~,v) - 0. The costs 
of these edges are defined as follows: c{vj ,v~j^-^) = p^{Sj) for 1 < j < n~. For each j we add the edge 
{■jsj 1 vj) with length and cost oo (for computational purpose a sufficiently large number M would do); this 
connects the node corresponding to the edge ej to vj that corresponds to Sj. See Figjsj 

The construction of T+ is quite similar except that the edge directions are reversed; assuming that the 
edges in S~^{v) are ordered such that £{ei, v) > £{e2,v) > . . . > £{e^+ , v), we create a path v — )• — )• 

. . ^ Vi with edge lengths £{e^+,v) — 0, ... , £{ej,v) — £{ej+i,v), . . . , £{ei,v) — £{e2, v). The costs for 
the edges in this path are set to pt{S^+), . . . , pt{Si) where Sj = {ei, . . . , ej}. For each j we add an edge 
(f^, 7ej ) with length and cost 00. This finishes the description of H. We now describe various properties 
of the graph H. Several of these properties are staright forward from the description of the construction and 
we omit proofs of the easy claims. 

The proposition below asserts the cost of the fractional solution in the edge-capacitated network H is the 
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Figure 3: Illustration of the reduction. Only Sq{u) and Sq{v) shown. The costs on edges in H are shown 
but not their lengths. The lengths of the infinite cost edges is and £'{u2,uf) = £{a,u) — £{e,u) and 
e{v^,v^)=i{c,v)-£{d,v). 

same as the cost of the fractional solution in the polymatroidal network G. 
Proposition 2. Eeei?« " ^'(e) = E.ey(r «) + P+(d+)). 

Proposition 3. For any edge e £ 5q (v) the length of the unique path in from the node 7e to v is equal to 
£{e, v). Similarly for e G 5q{v), the length of the unique path in from the node v to the node 7e is equal 
to £{e, v). 

We now establish a correspondence between paths in G and H that connect nodes in V. Let e = (n, v) 
be an edge in G. We obtain a canonical path q{u, v) from u to v in as follows: concatenate the unique path 
from u to 7e in with the unique path from 7e to v in T~. For any two nodes s,t € V let Vg{s, t) be the set 
of (simple) s-t paths on G and similarly Vh{s, t) be the paths in H. We create a map g : Vg{s, t) — )• Vh{s, t) 
as follows. Consider a path p G Vcis, t); we obtain a path p' G Vh{s, t) corresponding to p as follows. We 
replace each edge {u, v) e phy the canonical path q{u, v). 
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Lemma 4. The map g is a bijection. Moreover, for any two nodes u,v ^V, distii {u, v) = dist^^u, v). 

Now we establish a correspondence between cuts in G and H. For a given set of edges F C E' let 
sepQ(F) be set of node pairs inV x V separated by F in the graph G. Similarly for a set of edges F' C Eh 
let sepjj{F') be the set of node pairs inV x V separated by F' in the graph H. We say that a set of edges F 
is minimal with respect to separating node pairs if there is no proper subset of F that separates the same node 
pairs as F. 

Proposition 4. Let F' C Eh be minimal with respect to separating node pairs inVxV and of finite cost. 
Then for any v & V, F' contains at most one edge from and at most one edge from T^. 

Proof. Consider a node v and edge-sets F' n and F' n T+. For an edge e £ E there is a node G Vh 
and there is exactly one edge coming into 7e and exactly one edge going out of 7e and both are of infinite 
cost. Therefore, if F' is of finite cost, F' n consists of some edges in the path f j~ — )• . . . — )• v~_ — )• v 

fly 

contained in . Since the only way to reach v is through T~ it follows that if F' contains an edge {vj , 
then it is redundant to remove an edge {v^ , v~^^i) for i < j. Thus minimality of F' implies F' contains exactly 
one edge from . The reasoning for T+ is similar. □ 

Lemma 5. Let F' C Eh be minimal with respect to separating node pairs inV xV and of finite cost. There 
exists a set of edges F Q E such that sepQ{F) 5 sepniF') and ^{F) < c{F'). 

Proof. Given a minimal F' we obtain a set of edges F Q E as follows. From the proof of Proposition |4] we 
see that for any node v, F' contains at most one edge from and in particular if it contains an edge then it is 
an edge {v~ ,v~,-^) for some 1 < j < (for simplicity we identify v with ). Suppose there is such an 

edge e' = {vj , Vj^-^) in F' . Note that e' corresponds to the set — {ei, . . . , ej} of edges in ^q{v) ordered 
in increasing order by £{e, v) values. We add Sj to F and assign these edges to v in upper bounding z^(-F): by 
construction c(e') = p~{Sj). We do a similar procedure if e' E F n T+. It follows that the edge set F that 
we construct satisfies the property that v{F) < c{F'). 

We now show that sepg.(F) 5 sep^(F') Consider a pair (s, t) such that s is separated from t by F' in 
H. Suppose (s, t) is not separated by F in G. Let p be an s-t path that remains 'mG\F. From Proposition [5] 
there is a unique path g{p) G Vnis, t). For every edge e = {u,v) £ p consider the canonical path q{u, v) in 
H. Since e is not in F it implies that u can reach 7e in i/ \ F' and that 7e can reach v in H\ F'. This means 
that q{u, v) exists in if \ F'. This would imply that g{p) exists in H \ F' contradicting that assumption that 
(s, t) is separated by F'. □ 

We summarize the properties of the reduction. We assume that we have a polymatroidal network G = 
{V, E) with k demand pairs {si,ti), . . . , {sk, tfc) with associated demand values Di, . . . , D^. For all the 
cut problems of interest, the relaxations in Section |3] produce a length function i : E ^ and for each 
e = (n, v) associated non-negative values i{e, u) and £(e, v) such that l{e) = l{e, u) + ^(e, v). As before we 
use d~ and to denote the vector of ^(e, v) values for the incoming and outgoing edges at v. The reduction 
produces an edge-capacitated network H = {Vh, Eh) with the following properties: 

• each node of y is a node in Vh 

• for all u,v ^V, dist£(n, v) = dist^/ (u, v) 

• for any set of edges F' C Eh there is a corresponding set F C F such that sepQ(F) 5 sepji^(F') and 
u{F) < c{F'). 
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We also note that the reduction can be carried out in polynomial time. Moreover, given a set F' C Eh a 
set -F C that satisfies the last property in the list above can be found in polynomial time. 

We build on the reduction to obtain flow-cut gap results, all of which are based on using the relaxations 
from Section |3] which are dual to the corresponding flow problems. We argue via the reduction and known 
results on edge-capacitated networks that there exist integral cuts within some factor a of the fractional solu- 
tion. 

4.2 Multicut 

We consider the multicut problem for arbitrary demand pairs as well as symmetric demands. The relaxation 
satisfies the constraint that dist£(sj, tj) > 1 for each demand pair (sj, tj). The reduction from the preceding 
section produces a graph H = {Vh,Eh) and a fractional solution i' : Eh — M+ such that dist£/(sj, ti) > 1. 
We note that £' is a feasible solution for the standard distance based relaxation for multicut in edge-capacitated 
networks which is the dual for the maximum throughput multicommodity flow problem. The integrality gap 
of this relaxation has been studied and several results are known. Let /? = YIsgEh ^(^)^'(^) fractional 
solution value. Then one can obtain an integral multicut F' with cost c{F') that can be bounded in terms of 
/3. We summarize the known results. 

• Cheriyan, Karloff and Rabani ||8| showed that there exists an F' such that c(F') < 0(1) • /3^; this was 
improved by Gupta |[T6l to show the existence of a multicut F' such that c(F') < 0(1) • These 
results hold under the assumption that c(e) > 1 for all e. 

• Agrawal, Alon and Charikar ||2| improving the results in |[8l[T6l showed the existence of a cut F' such 
that c(F') = 0(niV23) . ^_ Here n is the number of nodes in the graph. 

• Saks, Samorodnitsky and Zosin 041 showed that there exist instances on which every integral multicut 
has a value Q{k) ■ /3. 

• Chuzhoy and Khanna |[9l showed that there exist instances on which every multicut has a value '') • 
/3. Further, they showed that the multicut problem is hard to approximate to within a factor of r2(2^°s " ") 
unless NP C ZPP. 

Since polymatroidal networks generalize edge-capacitated networks it follows that all the lower bounds in 
the above hold for the polymatroidal network case as well. The reduction also allows us to obtain upper bound 
for polymatroidal networks. We have to careful when using bounds that depend on the number of nodes in 
the graph. The reduction takes G with n nodes and m edges and produces an edge-capacitated graph H with 
n + 2m nodes. In the worst case H has il(n^) nodes. We thus obtain the following theorem. 

Theorem 1. In a directed polymatroidal network G on n nodes, for any given multicommodity flow instance 
with k pairs, if (3 is the maximum throughput multicommodity flow then: 

• There is a feasible multicut F' such that i^{F') < 0(1) • assuming that and p~ are integer valued 
for all V €z V. 

• There is a feasible multicut F' such that i'{F') < 0(n^^/^'^) • (3. 

Moreover, there exist polynomial-time algorithms to find multicuts guaranteed as above. 
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Symmetric demands: We now consider the symmetric demand case when a muhicut corresponds to sepa- 
rating {si,ti) or {ti, Si) for a given demand pair Sjtj. The relaxation for this has a constraint that dist£{si,ti) + 
dist^(tj, Sj) > 1. In contrast to the strong negative results for the general multicut problem, poly-logarithmic 
upper bounds on flow-cut gaps are known for symmetric demands in standard networks. In particular Klein 
et al. [23 1 show that if (3 is the cost of a fractional solution then there exists an integral multicut of cost 
0(log^ k) ■ (3. Even et al. [11] showed the existence of a multicut of cost O(lognloglogn) • /3. Note that 
these bounds are incomparable in that depending on the relationship between k and n one is better than the 
other. It is also known that there exist instances on which the gap is at least r2(log n). Via the reduction we 
obtain the following. 

Theorem 2. In a directed polymatroidal network G on n nodes, for any given multicommodity flow instance 
with symmetric demands on k pairs, the minimum multicut is 0(min{log^ fc, log n log log n}) • /3 where {3 is 
maximum throughput multicommodity flow for the symmetric demands. 

Remark 1. The flow-cut gap in polymatroidal networks for multiterminal flow^can be shown to be 2 via the 
reduction and the result ofNaor and Zosin / I?0|/ . 

4.3 Sparsest cut 

Now we consider the sparsest cut problem where the goal is to find a set of edges F to minimize iy{F) /D{F) 
where D{F) is the total demand of the pairs separated by F. The relaxation corresponds to finding edge 
length variables I to minimize the fractional cost subject to the constraint that Di ■ disti{si, ti) = 1. Via 
the reduction we produce an edge-capacitated network H such that Di ■ dist^/(sj, fj) = 1 and with the 
fractional cost preserved. In edge-capacitated networks there is a generic strategy that translates the flow-cut 
gap for multicut into a flow-cut gap for sparsest cut at an additional loss of an 0(log^^ Di) factor due to 
Kahale lfT9l (see also ll37l ): this has been refined via a more intricate analysis in 1311 to lose only an 0(log k) 
factor although one needs to apply it carefully. In ||2l a simple reduction that loses an O(logn) factor is 
given (this builds on |T9l). For directed graphs the known-gaps for sparsest cut are essentially based on using 
the corresponding gap for multicut and translating via the above mentioned schemes. We thus obtain the 
following results. 

Theorem 3. In a directed polymatroidal network G on n nodes, for any given multicommodity flow instance 
with k pairs, if (3 is the value of the maximum concurrent flow then there is a cut ofsparsity at most 0(n^^/^'^) • 
/3. 

Theorem 4. In a directed polymatroidal network G on n nodes, for any given multicommodity flow instance 
with symmetric demands on k pairs, there is a cut ofsparsity 0(min{log'^ k, log^ n log log n}) ■ f3 where (3 is 
maximum concurrent flow. 

5 Flow- Cut Gaps in Undirected Polymatroidal Networks 

In this section we consider flow-cut gaps in undirected polymatroidal networks. As we already noted, node- 
capacitated flows are a special case of polymatroidal flows. We show that line embeddings with low average 
distortion introduced by Matousek and Rabinovich [29| (and further studied in [32|) are useful for bounding 
the gap between the maximum concurrent flow and sparsest cut; we are inspired to make this connection from 
|[T3l who considered node-capacitated flows. For multicut we show that the region growing technique from 

■*In multiterminal flows we have a set of k terminals {si, S2, . . . , Sk} and flow can be sent between any pair of terminals; the goal 
is to maximize the total flow. The corresponding cut is referred to as multiterminal cut or multiway cut in which the goal is to remove 
a minimum-cost set of edges to disconnect every (ordered) pair of terminals. 
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ll25l that was used in lfT4ll for edge-capacitated multicut can be adapted to the polymatroidal setting. These 
techniques are also applicable to directed graphs — we defer a more detailed discussion. 



5.1 Maximum Concurrent Flow and Sparsest Cut 

We start with the definition of line embeddings and average distortion. 

Let {V, d) be a finite metric space. A map 17 : F — >• M is an embedding of V into a line; it is a contraction 
(also called 1-Lipschitz) if for all u,v ^ V, 

\g{u) -g{v)\ < d{u,v). 

Given a demand function w -.V xV ^ IR+ and a contraction : y — )• M, its average distortion with respect 
to w is defined as 

^ , , Y.u,v&vw{u,v)-d{u,v) 

2lu,v(^vw{u,v) ■ \g{u)-g{v)\ 
The following theorem is implicit in Q ; see lITSll for a sketch. 

Theorem 5 (Bourgain Q). For every n-point metric space [V^ d) and every weight function w :V y<V ^ M+ 
there is a polynomial-time computable contraction : y — t- M such that avgd^(g) = O(logn). Moreover, if 
the support ofw is k there is a map g such that avgd^{g) = 0(log k). 

Using the above we prove the following. 

Theorem 6. In undirected polymatroidal networks, for any given multicommodity flow instance with k pairs, 
the ratio between the value of the sparsest cut and the value of the maximum concurrent flow is 0{[ogk). 
Moreover, there is an efficient algorithm to compute an 0(log h) approximation to the sparsest cut problem. 



Recall the relaxation for the sparsest cut from Section |3.3| and the associated notation. To prove the 
theorem we consider an optimum solution to the relaxation and show the existence of a cut whose sparsity is 
0(log k) times the value of the relaxation. Let {V, d) be the metric induced on V by shortest path distances 
in the graph with edge lengths given hy I : E ^ M_|_ from the optimum fractional solution. Let g : V ^ M. 
be line embedding guaranteed by Theorem [5] with respect to d and the weight function given by the demands 
Di, that is w{si,ti) = Di for a demand pair and is for any pair of nodes that do not correspond to a demand. 
Without loss of generality we can assume that g maps V to the interval [0, /3] for some /3 > 0. For 6 G (0, /?) 
let Sg = {u \ g{u) < 9}. We show that there is a such that ^(S'e) is an approximately good sparse cut. Let 
D{6{S0)) be the total demand of pairs separated by Sg, that is D{6{Se)) = J2rSg separates Siti ^i- 



Lemma 6. 

D{d{Se))de = n 



1 

log k 



Proof. From the definition of D {6 (Sg)), 

/ D{6{Se))de= / ( V D,)d9 = y^Di - / lses.^.r.t..s^ud9 = y^Di-\g{si) - g{ti)\. 

Jo -/o • o . i -1^0 ■ 1 

f.oo separates Siti i=l i=i 

From the properties of g, 

A • \9[Si) - g{ti)\ 

We have the constraint J2i ^« " d{si,ti) = 1 from the LP relaxation; this combined with the above inequality 
proves the lemma. □ 
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The main insight in the proof is the following lemma. A version of the lemma also holds for directed 
graphs that we address in a remark following the proof. 



Lemma 7. 

u{5iSe))de<2y^Pu{du). 

Proof. Consider an edge uv G S{Se) and for simplicity assume g{u) < g{v). The length of e in the em- 
bedding is i'{e) = \g{v) — g{u)\ < l{e). The edge {u,v) G ^{Se) iff 9 is in the interval [g{u), g{v)]. 
Note that the cost v{5{So)) is in general a complicated function to evaluate. We upper bound v{5{Se)) 
by giving an explicit way to assign e = uv to either u or w as follows. Recall that in the relaxation 
l{e) = £{e, u) + i{e, v) where i{e, u) and £{e, v) are the contributions of u and v to e. Let r = ^jj^ and let 
i'{e, u) = r£'{e) and £'{e, v) = {1 — r)i'{e). We partition the interval [g{u), g{v)] into [g{u),g{u) + £'{e, u)) 
and [g{u) + £'{e, u),g{v)]; if 6 lies in the former interval we assign e to u, otherwise we assign e to v. This 
assignment procedures describes a way to upper bound u^d^Sg)) for each 9. Now we consider the quantity 
Jq v{5{Se))d9 and upper bound it as follows. 

Consider a node u and let Lu = {uv G 5{u) \ g{v) < g{u)} be the set of edges uv that go from u to the 
left of u in the embedding g. Similarly Ru = {uv G 5{u) \ g{v) > g{u)}. Note that and partition 
d{u). Let be the vector of dimension \S{u)\ consisting of the values £'{e, u) for e G 6{u). We obtain 
from by setting the values for e G -R„ to and similarly d^ from d'^ by setting the values for e G to 
0. Since < £'{e,u) < £{e,u) for each e G 5{u) we see that d'„ < du and (component wise) and hence 
dy < du and d^ < d„. Since pu is monotone we have that pu{d^) < Pu{<i-u) and /)„(d^) < Pu{du) (see 
Proposition [TJ. 

We claim that 

' iyi6{Se))d9 <Y,iPu{d{:) + pu{d^)), 

which would prove the lemma. 

To see the claim consider some fixed 9 and v'{6{Sg)). Fix a node u and consider the edges in 6{u) n Sg 
assigned to u by the procedure we described above; call this set Ag ^. First assume that 9 < g{u). Then the 
edges assigned to u by the procedure, denoted by Ag^u = {e ^ Lu \ 9 > g{u) — £'{e,u)}. Similarly, if 
9 > g{u), Ag^u = {e £ Lu \ 9 < g{u) + £'{e, u)}. From these definitions we have 

/3 



For a fixed node u, 



f v{6{Sg))d9 [ Pu{Ag,a)d9. 

UGV-^O 



I Pu{Ag^u)d9 = / Pu{Ag^u)d9+ / Pu{Ag^u)d9 
Jo Jq Jq(u) 



Let Lu = {ei, 62, . • . , Ch} where < £'{ei,u) < £'{e2,u) < . . . < £'{eh, u). Then 

/ Pu{Ag^u)d9 = \{£'{ej,u) - /(ej_i, u))p({ei, 62, . . • ,6^}) 
•^0 ,=1 



The right hand side of the above, is by construction and the definition of the Lovasz extension, equal to 



Pu{d^). Similarly, /f. . Pu{Ag^u)d9 = Pn(d^). □ 
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Remark 2. An examination of the proof of the above lemma explains the factor of 2 on the right hand side; the 
edges in 6{v) can be both to the left and right ofv in the line embedding and each side contributes pu{(i^) to 
the cost. This is related to the technical issue about undirected polymatroid networks where the flow through v 
takes up capacity on two edges incident to v. For directed graphs one can prove a statement of the form below 
where 6^{Sg) is set of edges leaving Sg. Notice that there is no factor of 2 since one treats the incoming and 
outgoing edges separately. 




u{5+iSg))de <Y,ipZ{dZ) + ptidt))- 



The above statement gives an embedding proof of the maxflow-mincut theorem for single-commodity directed 
polymatroidal networks and has other applications. 

We now finish the proof of Theorem |6] via the preceding two lemmas. 

. u{5{Se)) ^ fiu{S{Se))d9 

iTLin ^ — 

8mP)Di5iSe)) - j^D{5{Se))de 

< 2'Y^pu{du)-0{\ogk) = 0{\ogk)^p 

u u 

The above shows that the sparsity of Sq for some 9 is at most 0(log k) times which is the value 

of the relaxation. Given a line embedding g there are only n — 1 distinct cuts of interest and one can try all of 
them to find the one with the smallest sparsity. The efficiency of the algorithm therefore rests on the efficiency 
of the algorithm to solve the fractional relaxation, and the algorithm to find a line embedding guaranteed by 
Theorem [5| both have polynomial time algorithms and thus one can find an 0(logA;) approximation to the 
sparsest cut in polynomial time. 

Remark 3. Node-weighted flows and cuts/separators can be cast as special cases of flows and cuts in poly- 
matroid networks. Our algorithm produces edge-cuts from line embeddings in a simple way even for node- 
weighted problems — the v cost of the edge-cut automatically translates into an appropriate node-weighted 
cut. In contrast, the algorithm in [13] has to solve several instances of s-t separator problems in auxiliary 
graphs obtained from the line embedding. 

5.1.1 Sparsest Bi-partition Cut 

We worked with general edge cuts so far, but for certain applications, it is necessary to work with a special 
type of edge cut called a bi-partition cut. In an undirected polymatroidal network, an edge-cut F is said to be 
a bi-partition cut if there exists a set 5 C y such that F = 6g{S). In the case of edge-capacitated undirected 
networks, it is well known that for any multicommodity flow instance, there always exists a sparsest cut that 
is a bi-partition cut. This does not hold for polymatroidal networks, however, a factor 2 gap can indeed be 
shown between the sparsest cut and the sparsest cut restricted to bi-partition cuts; moreover this factor is tight. 

Theorem 7. Given any edge cut F for a multicommodity flow instance in an undirected polymatroidal net- 
work G = {V,E), there exists a bi-partition cut 6g{S) whose sparsity is atmost 2 times the sparsity of F. 
Furthermore this factor is tight. 

The proof of the above theorem can be found in Section [C] of the appendix. Theorem [6] and Theorem. [7] 
together imply a logarithmic gap between maximum concurrent flow and sparsest bi-partition cut. This is 
formally stated in the following corollary. 

Corollary 1. In undirected polymatroidal networks, for any given multicommodity flow instance with k pairs, 
the ratio between the value of the sparsest bi-partition cut and the value of the maximum concurrent flow is 
0(log/c). 
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5.2 Maximum Throughput Flow and Multicut 

We prove the following theorem in this section. 

Theorem 8. In undirected polymatroidal networks, for any given multicommodity flow instance with k pairs, 
the ratio between the value of the minimum multicut and the value of the maximum throughput flow is 
Oi}.ogk). Moreover, there is an efficient algorithm to compute an Oilogk) approximation to the minimum 
multicut problem. 

We recall the relaxation for the minimum mulitcut problem from Section |3.2[ Consider an optimum 
solution to the relaxation given by edge lengths l{e), e ^ E and the partition of i{e) for each e = uv between 
u and V given by the variables £{e, u) and i{e, v). We will show that there exists a multicut F C E' for the 
given pairs such that = 0(log k){Y^^ pv{dy)). 

By slightly generalizing the proof of Lemma |7] we obtain the following. 

Lemma 8. Let g : V ^ [0, (3] be a contraction, let < ao < a < b < bo < 13 and Se = {u \ g{u) < 9}. 
Suppose for every edge e = uv £ Uo^^afil^iSe), g{u) and g{v) are both in [oq, 6o]- Then, 

J v{5{Se))de<2 Yl P-^^^-y 

° v:g{v)£[ao,bo] 

Proof. The proof is very similar to the proof of Lemma[7| except that to upper bound the left hand side in the 
statement of the lemma, we only need to consider edges that are in the set U0^[a.b]^{Se)- The condition in the 
lemma assures us that any node that is involved in 6{S0) have to lie within the interval [oq, bo]. Thus, it is 
sufficient to consider the set of nodes v : g{v) € [oq, 6o] in the integral on the right hand side. The proof is 
written out in detail in Sec. |B1 □ 

Given a graph G with edge lengths ]R_|_, a node v and radius r, let Bq{v, r) = {u\ dist£(t>, u) < 

r} denote the ball of radius r around v according to edge lengths £. We omit £ and G if they are clear from 
the context. For a set of nodes X C V we let vol(X) = J2vex Pvi^v) denote the total contribution of the 
nodes in X to the objective function. 

Lemma 9. Let 6 < 1 and suppose i{e) < 2\ogk '^^^ ^- Then, for any given node s and k > 2 there exists 
ar £ [0, 6) such that u{6{B{s,r)) <alogk- j(vol(i?(s, r)) + vol(y)/k), with a = 28. 

Proof. For simplicity we assume here that log k is an integer multiple of 3. Order the nodes in increasing 
order of distance from s: this produces a line embedding gg : V ^ M+. For integer i > define = 2^^- 
Define ao = vol(y)/A; and for i > 1 let ai = ao + wol{B{s, Vi)). 

Consider any 1 < j < 2 log k. We apply Lemma[8]to the embedding gs and the interval [rj_i, r^]; note 
that i{e) < which implies that we can indeed apply the lemma. Also any edge e G '<->ee[rj^i,rj]^{Se) 

satisfies the property that g{u) G [rj-2, rj+i] and g{v) G [rj-2, rj+i] since i{e) < Y^gk- Thus 

i^{6{B{s,e))de < 2 Yl Mdv) 

< 2(Qj+i - aj-2) (5) 

We claim that there is some 1 < j < 2 log k such that aj+i < 8aj_2- Suppose not, then 03, > 8a3(j_i) 
for all 1 < i < This implies that a^i > 8'ao = 2^'ao. Therefore, with i = this implies that 

a2iogfc > 22i°g^^^ > 4vol(y) which is impossible. 



/ 
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Thus there exists a j such that Oj+i < 8aj_2- Consider that j, equation ([5]) impUes that 



v{5{B{s,e))de < 2{aj+i-aj-2) 
< 2(7a,--2). 



If we pick r uniformly at random from the interval [rj_i,rj], where satisfies the above property, the 
expected cost of v{5{B{s, r))) is 

v{5{B{s,e))de < 



J-1 



from the preceding inequality and the fact that rj — rj^i = jt^a-- Hence there exists an r € [rj_i, r^] such 
that iy{6{B{s, r))) < ^^^aj^2- Since aj^2 - ao < vol(S(s, r)), the lemma follows. □ 

Now we consider the following algorithm for finding a multicut from a given fractional solution. 

• LetF^{eK(e)> j^i^}. 

• G' ^ G[E\F]. 

• While (there exists a pair Sjtj connected in G') do 

- Let Sjtj be a pair connected in G'. 

- ViaLemma|9]with5 = l/2findr < 1/2 such that i/((5G"(5G'(si) ^))) < 2alog A;-(vo1(Sg'/(sj, r))+ 



vol(y)/A;). 
_ F ^ Fu6G'{BG'{sj,r)). 

- Remove the vertices Bq' {sj,r) and edges incident to them from G'. 

• Output F as the multicut. 

Lemma 10. The set of edges F output by the algorithm is a feasible multicut for the given instance. 

Proof. (Sketch) One can prove this by induction on the number of steps in the while loop. We consider the 
first step. The diameter of the ball BQt{sj,r) is 2r < 1 and hence the end points of any pair cannot both 
be inside this ball. We remove the edges 6{BG'{sj,r)) and by the preceding observation there is no need 
to recurse on this ball. The algorithm recurses on the remaining graph G' — Bc'isj, r), and by induction 
separates any pair with both end points in that graph. □ 

Now we argue about the cost of the set F output by the algorithm. Let Fq <— {e \ £{e) > jj^^} be the 
initial set of edges added to F and let Fj be the set of edges added in the i'th iteration of the while loop. 

Lemma IL u{Fo) < 81og/c • I]^, /0^,(d^,). 

Proof. For v e V let Ay = {e e S{v) n Fq | £{e, v) > g-j^}- We can upper bound iy{Fo) by Pvi^v) 
since the latter term counts each edge uv e Fq in at least one of Au and A^ since £{e, u) + l{e, v) = £{e) > 
jj^p; . From the definition of the Lovasz extension 

.1 /-l/CSlogfc) ^ 

p,(d,) = / /9,(d^)d0 > / Pv{d'^y)de > - — -p^iAy), 

Jo Jo 8 log k 

where we used non-negativity of py for the first inequality above and monotonicity for the second. □ 
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Lemma 12. J2^>l ^i^i) < 4a log k 

Proof. (Sketch) From the algorithm description, Fj = 6{BG'{sj,r)) for someterminal Sj andradius r < 1/2 
where G' is the remaining graph in iteration i. Moreover, z^(Fj) < 2a log k • {vo\{BQ/{sj,r)) + vol(y)/A;). 
Since the nodes in Bqi {sj,r) are removed from the graph, a node u is charged only once inside a ball. Hence 

Y.'^i^i) < J^2alogfc- vol(y)/fc + 2alogA;^p^(d„) < 4a log A; ^ p^d^), 

i i V V 

since there are at most k iterations of the while loop; each iteration separates at least one pair. □ 
Since u is subadditive (see Lemma [T) 

iy{F) < i/(Fo) + ^zv(Fi) < (8 + 4a)logA:^/54d^). 

i>l V 

This finishes the proof of Theorem [8] 

6 Conclusions 

We considered multicommodity flows and cuts in polymatroidal networks and derived flow-cut gap results in 
several settings. These results generalize some existing results for the well-studied edge and node-capacitated 
networks. We briefly mention two results that can be obtained via the line embeddings technique that we did 
not include in this paper. A multicommodity flow instance in an undirected network G = {V,E) is a product 
multicommodity flow instance if there there is a non-negative weight function tt : V ^ M+ and the demand 
Duv between u and v is 7t{u) ■ 7r(v). The associated cut problem is interesting because it corresponds to 
finding sparse separators in graphs which in turn can be used to find balanced separators; these have several 
applications. It was shown in ||22ll that in edge-capacitated undirected planar networks, the flow-cut gap for 
product multicommodity flow instances is 0(1) (in fact they showed this holds for any class of graphs that 
excludes a fixed graph as a minor). Rabinovich f32l showed that the main technical theorem in [22] also 
leads to a line embedding theorem, and this was used in |13| to show an 0(1) flow-cut gap for product 
multicommodity flow instances in noJe-capacitated planar graphs. Our work here shows that this is true for 
undirected planar polymatroidal networks. Arora, Rao and Vazirani Q gave an O(vTogn) -approximation, 
via a semi-definite programming relaxation, for the sparsest cut problem in an undirected edge-capacitated 
network. Note that this is not a traditional flow-cut gap result since the SDP-based relaxation used is strictly 
stronger than the dual of the multicommodity flow relaxation. By interpreting the main technical result in 111 
as a line-embedding theorem, |[T3l obtained an O(vTogri) -approximation for sparsest cut in node-capacitated 
graphs; this can also be extended to the polymatroidal setting. 

Flow-cut gap questions for node-capacitated problems are less well-understood than the corresponding 
questions for edge-capacitated problems; line-embeddings provide a tool to obtain upper bounds on the gap 
but they do not provide a tight characterization as -embeddings do for the edge-capacitated case. We hope 
that polymatroidal networks and their applications to network information flow provide a new impetus for 
understanding these questions. 
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A Proof of Lemma [13] 



Lemma 13. For a directed polymatroidal network, the dual of the maximum throughput flow problem is 
equivalent (in terms of value) to the program given by Q (for undirected the program given by Qj. 

Proof. We will show the proof for the undirected case, the proof for the directed case is similar. The program 
for maximum throughput flow is given by: 

max ^ ^ f{p) 

S.t. 

e-.e&S P-e£p 

f{p) > VpeP(,^,t^),Vi = i...fc. 

The dual of the flow linear program can now be written. Let the dual variables dy{Sy) correspond to the 
non-trivial constraint in the above linear program. Then the dual linear program is: 



7'rf:=min ^ ^ dy{S)py{S) 
vev sc5{v) 

s.t. 

Yl ^«(^)+ E ^-(^) ^ 1 Vp e where e 

e=uv:e<^p \SCSiu):eeS SC5(t)):eeS' / 

du{S) > OVn G y V5 C 

This can be rewritten equivalently as 

Vd := min ^ ^ dy{S)py{S) 

v&V SCS{v) 

S.t. 

\SC5{u):eeS SCS{v):eeS 

dist£{si,ti) > 1 I <i < k 



uv 
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Let us define new variables i{e, u), £{e, v) for eacli edge e = uv, and rewrite the linear program: 



vev sc5{v) 

s.t. 

£{e) := £{e,u) + i{e,v), where e = uv 
£{e,u) = Y2 du{S) \/e £ E, e = uv 

SCS(u):eeS 

£{e,v) = Y2 dy{S) \/e£E,e = uv 

SCS{v):e£S 

dist£(sj,tj) > 1 l<i<k 
du{S) > 
£{e,u),e{e,v) > ^u e V \/S C S{u). 

The minimization is over the variables ^(e,n) a.nddy{S). Observe for any fixed t; the variables 5), S" C 
(5(f) influence only the variable £{e,v),e G S{v). Hence, for any v and a fixed assignment set of values 
£{e, v),e G d{v) the optimal choice of variables d,o{S), S S{v) can be obtained by solving the following 
linear program: 

mill ^ dy{S)py{S) 

SCS{v) 

s.t. 

Y2 dv{S) = £{e,v) Ve G e = 

SC5(v):e£S 

du{S) > 0, 5 C 6{v), Vf G V. 

Recalling the definition of the convex closure of a function, one sees that the value of the above linear program 
is equal to p^,(d^); note that for polymatroids we can drop the constraint dv{S) = 1 in the linear program 
for the convex closure. Since the convex closure is equval to the Lovasz extension we obtain the desired 
equivalence of the formulations. □ 



B Proof of Lemma |8] 

Lemma 14. Let g : V ^ [0, f3] be a contraction, let < ao < a < b < bo < P and Sg = {u \ g{u) < 9}. 
Suppose for every edge e = uv £ ^e£[a,b]^iSe)> di^) and g{v) are both in [oq, 6o]- Then, 

j\[5{Se))de<2 Yl P-i^^-y 

" t):g(t))e[ao,feo] 

Proof. Consider an edge uv G 6{Sg) and for simplicity assume g{u) < g{v). The length of e in the embed- 
ding is £'{e) = \g{v) — g{u)\ < £{e). The edge (n, v) G ^{Se) iff 9 is in the interval [g{u), g{v)]. Also by the 
conditions of the theory for every such (u, v), g{u) G [oq, 6o] and g{v) G [ao, 6o]. Note that the cost v{5{Se)) 
is in general a complicated function to evaluate. We upper bound v{5{Sq)) by giving an explicit way to assign 
e = uv to either u or w as follows. Recall that in the relaxation £{e) = £{e, u) + £{e, v) where £{e, u) and 
£{e, v) are the contributions of u and v to e. Let r = and let £'{e, u) = r£'{e) and £'{e, v) = {l — r)£'{e). 
We partition the interval [g{u), g{v)] into [g{u),g{u) + £'{e,u)) and [g{u) + £'{e,u), g{v)]; if 9 lies in the 
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former interval we assign e to u, otherwise we assign e to v. This assignment procedures describes a way to 
upper bound u{6{Sg)) for each 9. Now we consider the quantity J^^ i'{6{Sg))d9 and upper bound it as follows. 

Consider a node u and let Lu = {uv G 6{u) \ g{v) < g{u)} be the set of edges uv that go from u to the 
left of u in the embedding g. Similarly Ru = {uv G 5{u) \ g{v) > g{u)}. Note that Lu and Ru partition 
5{u). Let (^'u be the vector of dimension \5{u)\ consisting of the values £'(e, u) for e G (5(n). We obtain 
from d'jj by setting the values for e G i^^j to and similarly d^ from d'„ by setting the values for e G L„ to 
0. Since < £'{e,u) < i{e,u) for each e G 6{u) we see that d'„ < d^ and (component wise) and hence 
d^ < d„ and d^ < d„. Since pu is monotone we have that Pui^u) < Puidu) and ^^(d^) < ^^(dM) (see 
Proposition [TJ. 

We claim that ^ 

1^ v{6{Se))d9 < ^P-^'^u) + Pu{dZ))^ 

" tieV:3(w)e[ao,bo] 

which would prove the lemma. 

To see the claim consider some fixed 9 and v{5{Sq)). Fix a node u and consider the edges in 5{u) n Sq 
assigned to u by the procedure we described above; call this set Aq^u- First assume that 9 < g{u). Then the 
edges assigned to u by the procedure, denoted by Aq ^ = {e £ Lu \ 9 > g{u) — i'{e,u)}. Similarly, if 
9 > g{u), Ag^u = {e ^ Lu \ 9 < g{u) + l'{e, u)}. From these definitions we have 

JieV:g(M)e[ao,6o] 

^ l\{5{Se))d9 < Y f PuiAe,u)d9 

< Y Pu{Ae,u)d9. 

u&V:g{u)&[ao,bo\ 



For a fixed node u, 



/3 rgiu) /•/? 

= / pu{Ae^u)d9+ / , 
Jo Jg{u) 



Let Lu = {ei, 62, • • . , e/i} where < £'{ei,u) < £'{e2,u) < . . . < i'{eh, u). Then 

/ = 2^(^'(ej,u) - i'{ej-i,u))p{{ei,e2,. . .,ej}) 



The right hand side of the above, is by construction and the definition of the Lovasz extension, equal to 

ig{u) ' 



Pu{d^). Similarly, /jj . Pu{Ae,u)d9 = /5„(d^). □ 



C Proof of Theorem HI 

We recall the statement of Theorem |71 

Theorem. Given any edge cut F for a multicommodity flow instance in an undirected polymatroidal net- 
work G = {y, E), there exists a bi-partition cut 6g{S) whose sparsity is atmost 2 times the sparsity of F. 
Furthermore this factor is tight. 
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Proof. Let Vi,V2, ■ ■ ■ ,VhhQ the connected components of G induced by the removal of the edge-cut F. Let 
D{F) be the total demand separated by F. We show a cut F' = 5g{S) such that D{F') > D{F)/2 and 
F' C F. 

We obtain F' as follows. Construct an undirected graph with nodes vi,...,Vh, corresponding to the sets 
Vi,. .. ,Vh. For each pair Vi, vj we add an edge ViVj with weight w{viVj) equal to the total demand of all 
pairs with one end point in Vi and the other end point in Vj. Note that the total weight of all edges is equal to 
D{F). It is well-known that in any undirected weighted graph there is a partition of the nodes into A and A'^ 
(the complement of A) such that the total weight of edges crossing the partition is at least half the weight of all 
the edges in the graph (a random partition guaratees this in expectation and gives a simple 1/2-approximation 
to the NP-Hard maximum-cut problem). Let A be such a cut in H; we have w{5h{A)) > D{F)/2. Now 
consider the set of nodes S = U^.g^Vi in G and the corresponding cut F' = 5g{S). It follows that D{F') > 
D{F)/2. Moreover F' C F and by monotonicity of i^{F') < v{F). Hence, 

D{F') - D{Fy 

which implies that the sparsity of F' is at most twice that of F. By construction F' is a vertex bi-partition cut. 

To see that the factor of 2 is tight, consider a polymatroidal network G = {V,E) which is a star on n 
nodes {vq, ui, . . . , Vn-i} with center vq; the edges are ei, . . . , e^-i where = VQVi. Assume n is even. 
The only capacity constraint is a polymatroidal constraint at node vq, which constrains the total capacity of 
every subset of {ei, e„_i} by a value of 1. The demand graph is a complete graph on the nodes with 
each pair having a unit demand. Now consider an edge cut F which removes all the edges: i^(-F) = 1 and 
D{F) = (2), so the sparsity is ^(^^_iy whereas any bi-partition cut F' = S{S) also has value i^{F') = 1 and 

D{F') = |S'||S'|'^, which means the sparsity is minimized with l^l = | and is given by The ratio of the 
two sparsities is 2(1 — ^) and approaches 2 as n — > 00. □ 
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